Giant Relaxation Oscillations in a Very Strongly Hysteretic SQUID ring- Tank Circuit 
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In this paper we show that the radio frequency (rf) dynamical characteristics of a very strongly 
hysteretic SQUID ring, coupled to an rf tank circuit resonator, display relaxation oscillations. We 
demonstrate that the the overall form of these characteristics, together with the relaxation oscil- 
lations, can be modelled accurately by solving the quasi-classical non-linear equations of motion 
for the system. We suggest that in these very strongly hysteretic regimes SQUID ring-resonator 
systems may find application in novel logic and memory devices. 

PACS numbers: 85.25.Dq, 74.50. + r 



INTRODUCTION 



Over the years the non-linear, dynamical properties 
of SQUID rings (i.e. thick superconducting rings con- 
taining either one or two Josephson weak link devices) 
have formed the basis of a range of technologically impor- 
tant devices [jj , most recently through the phenomenon 
of stochastic resonance 0. Apart fr om these applica- 
tions, the studyof the physics of the SQUID rings them- 
selves 0, 0, [J EI has added considerably to the general 
field of non-linear dynamics. As studied previously, these 
dynamical properties have been due to the quasi-classical 
behaviour of SQUID rings in the regime where quantum 
processes can be neglected. However, of late much in- 
terest has been shown in the use of SQUID rings as 
purely quantum devices in possible quantum technolo- 
gies This interest has been encouraged by recent ex- 
perimental work on superp osition states in weak link cir- 
cuits B^JSIISIIlllllIl and, in particular, on SQUID 
rings [ljjlla]- Inevitably, at some stage, quantum circuits 
are required to interact with classical probe circuits if in- 
formation about their quantum state/evolution is to be 
extracted. As we have shown, for the case of the SQUID 
ring this interaction leads to the growth of non-linear be- 
haviour in the classical part of the system 0, 0, 0, 0] . 
Thus, from this perspective, non-linear dynamics always 
has a crucial role to play in the description, quantum 
or quasi-classical, of SQUID rings. In this paper we shall 
demonstrate that, even with the quite understandable at- 
tention being paid to SQUID rings at the quantum level, 
there is much still to be explored in the quasi-classical 
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FIG. 1: Schematic of the coupled SQUID ring-rf tank cir- 
cuit system also showing the static magnetic flux bias Qxstat 
applied to the ring. 



regime. This can add greatly to our knowledge of non- 
linear systems and could lead to novel and powerful ap- 
plications. 

In the description of SQUID rings [jj it is customary 
(from this quasi-classical viewpoint) to identify regimes 
of behaviour in terms of the parameter f3 (= 27rA/ c /$o), 
where A s is the ring inductance, I c is the (Josephson) 
critical current of the weak link in the ring and $o = 
h/2e. In the presence of an external magnetic flux $> x the 
response of a SQUID ring is to set up a screening super- 
current 7 S ($ x ). The functional form of I s (& x ) depends 
on the value of (3. Thus, for (3 < 1, I s is always single 
valued in <fr x while for (3 > 1 this current is multi-valued 
and hysteretic in Q x . In this latter regime I s (Q x ) can be 
viewed as a set of approximately diamagnetic branches 
with switching between adjacent branches (in the ab- 
sence of noise) occurring when I s reaches ±7 C . Given this 
branch switching, hysteresis loops will be swept out, with 
concomitant dissipation, when I s exceeds ±/ c . Devices 
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which utilize the properties of single weak link SQUID 
rings, such as magnetometers pj, are usually operated at 
low [3 values (e.g. 2 to 5). In operation as magnetome- 
ters single weak link SQUID rings are typically coupled 
inductively to a radio frequency (rf w 20MHz), paral- 
lel LC, tank circuit which is excited using an rf current 
source. This arrangement is shown schematically in fig- 
ure 1. As is well known, the rf voltage developed across 
the tank circuit is a function of both the rf magnetic flux 
amplitude and the static magnetic flux (<fr xs tat) applied 
to the ring. This static flux dependence is periodic in $0, 
hence the use of the SQUID ring-tank circuit system in 
magnetometry. 

The dynamical behaviour of a SQUID ring coupled to 
an external resonant circuit is governed, ultimately, by 
the form of its potential U (<&, $ x )- Quasi-classically [lj 
(i.e. when the effective mass of the ring - the weak link 
capacitance C - is large, typically 10 -12 to 10 _14 F for 4 
Kelvin operation) this potential has the form 

where <E> (= <§> x + AI S ) is the total included flux in the 
ring. Clearly, if $> x changes slowly with time (i.e. slowly 
compared with any ring time constants), the ring will 
always remain at, or very close to, the minimum in 
U (<&, $ x )- In dealing with ring- resonator dynamics it has 
been customary to assume that this adiabatic constraint 
is satisfied, e.g. in the low (3 rf SQUID magnetometer 
where the peak rf flux at the ring typically > $0 (= (J-f 
for a fraction /1 of the tank circuit flux ip coupled to the 
ring) . The current I s set up in response to this flux then 
couples back to the tank circuit, and so on. This back re- 
action, which reflects I s ($ x )i affects the rf voltage V ou t of 



the tank circuit in a non-linear way. With Ii n linearly am- 
plitude modulated in time (e.g. at w 100Hz), the result is 
the rf SQUID magnetometer (V ou t verses h n ) character- 
istics. Following an initial linear riser in Vout versus Ii n , 
these characteristics consist of a set of voltage plateaux 
(steps) spaced at regular intervals along ij„, with pe- 
riod cc $o/A. The <3>o-periodic modulation of these step 
features by a static (or quasi-static) applied flux $> X stat 
forms the basis of the rf SQUID magnetometer. 

In this paper we report on the experimental and the- 
oretical investigation of single weak link SQUID rings 
in the so-termed very strongly hysteretic regime, a pre- 
viously little explored part of the parameter space of 
SQUID behaviour. With the SQUID ring inductively 
coupled to a high quality factor (Q) radio frequency 
(rf w 20MHz) resonant circuit, we show that quite re- 
markable non-linear phenomena (large amplitude relax- 
ation oscillations) can develop. We demonstrate that ex- 
actly these dynamics are solutions of the coupled non- 
linear equations of motion of the system. 
II. SQUID RING- TANK CIRCUIT DYNAMICS 

IN THE STRONGLY HYSTERETIC REGIME 

A. Quasi-classical equations of motion 

The quasi-classical description of the ring-tank circuit 
system is based on the resistive ly shunted junction plus 
capacitance (RSJ+C) model of a weak link. Here, the 
supercurrent channel through the link is in parallel with a 
normal current channel (resistance R s ) and a capacitance 
C s . In the model the normal channel turns on when the 
link current I > I c . With representative values of C s 
from 10~ 12 to 1CU 13 F, generally accepted resistances R s 
are in the 10 to 100O range (l9|. With this in mind, the 
coupled ring-tank circuit equations of motion are 0, |20| 



SQUID ring : C s ^ + ±§ + 7 c sin ( *g) + *_ = (2) 
dt 2 R s dt V *o / A s (1 - K 2 ) A s (1 - K 2 ) 

„ d 2 in 1 dip ip , . £t$ , . 

Tank circuit : C T J + -J + - r2) = I m (t) + jrf^ (3) 



for tank circuit capacitance, inductance and parallel res- 
onance resistance Ct, Lt and Rt, respectively, coupling 
constant K = yj M 2 /LtA. for mutual inductance M, and 
Josephson current I = I c sin (27r$/$o)- It is apparent in 
that there are several possible ring time constants 
(C S R S ,A S /Rs, 1/y/AgCg). Clearly, if the adiabatic con- 
straint holds, and there are no time constant problems, 
the derivative terms in @ can be neglected and © re- 
duces to 



J c sin (2tt$/$ ) + $/A s (1 - K 2 ) = ^/A s (l - K 2 ) 

(4) 

In the past this, together with J5J, has been the starting 
point of most dynamical descriptions of ring-tank circuit 
systems, at least for low fj SQUID rings. We note that 
in this linearized, small (3, description I c (i.e. (3) can be 
found from the point of onset in I in of the periodic step 
features in V ou t versus Ii n (the breakpoint) beyond the 
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initial linear riser |21| . 

As we emphasized recently |2jj, the adiabatic con- 
straint breaks down when the ring is strongly hysteretic 
and underdamped p^ . Qualitatively, at 4 Kelvin, transi- 
tion widths for branch switchings in I s ($> x ) are typically 
« $o/100. At 20MHz, and at the rf drive levels used ex- 
perimentally for a (3 ss 50 — 100 ring (fnp > 100<&o), the 
available time to switch between branches w 10~ 12 secs, 
comparable to the shortest ring time constant. In this 
range of (3 novel multiple level structures develop in the 
experimental V ou t versus U n characteristics |20j . These 
structures (plateaux of almost constant V out along I in ) 
are on a scale in Ii n large compared with standard rf 
SQUID magnetometer steps. However, these still main- 
tain a form of $o-periodicity in Q xs tat- We also demon- 
strated numerically that these multiple level structure 
characteristics formed a new and particular class of so- 
lutions of ( I2J) and ■ Given the range of non-linear 
phenomena that have already been recognized in ring- 
tank circuit systems, it would be surprising if these mul- 
tiple level structures constituted the only new class of 
solutions in the high (3 regime. At even higher (3 a sec- 
ond type of non-linear dynamics (also a solution of|21and 
© develops - relaxation oscillations in rf voltage against 
current - on scales very large (« x 10) compared with 
conventional SQUID magnetometer characteristics. 



B. Experimental V ou t versus 7j n characteristics 
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FIG. 2: Experimental (4.2 Kelvin, <I>xstat = n<E> , A = 6 x 
KT 10 H, L t = 6.3 x 10" 8 H, / If = 22.671MHz, K 2 = 0.008 
and Q — 515) K ou t versus I- m characteristic for a niobium point 
contact SQUID ring-rf tank circuit system in the very highly 
hysteretic regime showing large scale relaxation oscillations. 
Here, the triangular ramp amplitude modulation of the rf is 
at 18Hz and the detection bandwidth of the receiver is 1MHz. 



In our experiments we made use of Zimmerman, 2-hole, 
niobium point contact SQUID rings [23. The point con- 
tact weak links in these rin g's w ere adjusted, in situ, at 
liquid helium temperatures [Til l20| . a technique which 
allows us to make weak links in a controlled way over a 
very wide range of (3 values (e.g. 1 < (3 < 100 ). In our 
experimental arrangement the tank circuit voltage V ou t 
(figure 1) was first amplified by a liquid helium cooled 
GaAsFET preamplifier (gain ~ 20dB, noise temperature 
< 10 Kelvin). This signal was further boosted using a re- 
ceiver of our own design and then diode detected. Over- 
all, the receiver system (including the GaAsFET pream- 
plifier) was designed to combine very low noise with very 
large dynamic range and slew rate, all three properties 
being required to observe the very strongly hysteretic be- 
haviour reported here. 

In figure 2 we show a perfectly typical experimen- 
tal diode detected V out versus Ij„ SQUID magnetome- 
ter characteristic in this very strongly hysteretic regime. 
This was taken at 4.2 Kelvin, in a bandwidth of 1MHz, 
with Q = 515, K 2 = 0.008 and a static bias flux $> xs tat = 
n<&o, n integer, supplied by a second coil coupled to the 
ring (figure 1). Here, adopting common practice, this 
characteristic was recorded at the n$ Q bias frequency f r f 
(22.671MHz in this case) equal to the resonant frequency 
of the ring-tank circuit system when <fr xs tat = n&o- As 
is standard for these dynamical (SQUID magnetometer) 
characteristics, the rf is amplitude modulated using a low 



frequency triangular ramp. In this example the ramp fre- 
quency is 18Hz. As can be seen, the characteristic con- 
sists of a first feature, a flat step, followed by two further 
step features in which it is apparent that V ou t oscillates 
very rapidly as increases. Although not shown in fig- 
ure 2, these oscillatory step structures, periodic in Jj ra , 
continued to repeat as Jj n was made larger. The rf volt- 
age and current scales in figure 2 are made clear by the 
small, subsidiary characteristic in the same figure. This 
was taken at 4.2 Kelvin, with ^ xs tat — n^o, using the 
SQUID ring and tank circuit of the main characteris- 
tic but with the weak link in the ring adjusted to give 
a (3 just greater than unity, i.e. the regular rf SQUID 
magnetometer regime. It is obvious that the scalings of 
the main characteristic are very large indeed compared 
to those for a standard SQUID magnetometer. In fig- 
ure 3 we show details of the second step structure of 
figure 2 along which the large scale oscillations first ap- 
pear. Here, we have expanded the first section of this 
step structure (shown boxed in figure 3). In the expanded 
picture it is apparent that V ou t is undergoing rapid os- 
cillations as Ii n is increased. Furthermore, the spacing 
in Ii n between adjacent oscillations decreases monotoni- 
cally as U n is swept across the step; this continues right 
across the step (not shown in figure 3). We note that 
the effect of a bias flux Q xs tat on these oscillatory step 
structures was minimal, only changing the amplitude of 
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FIG. 3: Expanded region of the second plateau in figure 2 
showing details of the relaxation oscillations. 



the oscillations very weakly, with a periodicity of $o- Ex- 
perimentally, a search of the first (flat) step in figure 1 
did not reveal any voltage oscillations, only noise. Again, 
the effect of $ X stat was found to be negligible. 

The waveform of figure 3 reflects a limit cycle consist- 
ing of a very slow build up (at least on the time scale 
of one rf oscillation 27r/o; r /) followed by an extremely 
fast discharge with energy drawn from the tank circuit. 
This is the followed by the next slow build up with en- 
ergy supplied from the rf current source over many rf 
periods, and so on. These are, of course, relaxation oscil- 
lations which can be found in many ph ysic al situations 
where non- linear dynamics is involved [28L 12^. I30L l3l| . 
For the particular case of the driven SQUID ring-tank 
circuit system of figure 3, the peak rf amplitude is in- 
creasing (linearly ramped with time) from left to right. 
Since the energy is made up from the rf current drive, this 
means that the time interval between successive oscilla- 
tions shortens as ij n increases, as is apparent in figure 3. 

Without the combination of very low noise receiver 
(including a liquid helium cooled GaAsFET preamplifier 
stage) and very large dynamic range/slew rate capability, 
the relaxation oscillations shown in figures 2 and 3 would 
not be observed. For example, reduced bandwidth or 
added noise, or a combination of both, can wash out 
these oscillations. What then remain are large (compared 
with ordinary SQUID magnetometer characteristics - see 
insert figure 2), roughly constant rf voltage, current steps. 
This is presumably why these dynamics have not been 
observed in the past. 

It is interesting to note that the relaxation oscilla- 
tions observed by us in coupled SQUID ring-tank cir- 



cuit systems in the very strongly hysteretic regime ap- 
pear to have analogues in the phase slippage experi- 
ments which have been performed on superfluid helium 
systems [2jJ H(J |27|]. It seems clear that, putting the 
fundamental description of superfluid and superconduct- 
ing condensates aside, the phase/flux slippage in the two 
systems is analogous. At a phenomenological level this 
means that the equations of motion governing the dy- 
namics of the two ring-resonator systems are also analo- 
gous. In the highly hysteretic regime for phase slippage 
in superfluid rings containing weak link orifices we would 
therefore expect to see (giant) relaxation oscillation phe- 
nomena and, perhaps, multilevel dynamics |20|. 



C. Numerical calculations 

The Vout versus /;„ characteristics of figures 1 and 2 
can be modelled using and © which can handle the 
problem of non-adiabaticity perfectly well |20j. However, 
since for typical SQUID ring parameters (C s ps 1CU 13 F, 
A s = 6 x 1CU 10 H, R s ps 10fi), the ring time constants 
will be 2 or more orders of magnitude shorter than the 
tank circuit time constant (= 1 /20MHz) , great care must 
be taken to ensure that accurate solutions are found. To 
solve (J2J and © we used fourth order Runge-Kutta with 
an adaptive step size algorithm 0, 0, [2(J . We were also 
careful allow a sufficient number of steps per rf cycle to 
follow the changes in the SQUID ring accurately. For our 
comparison with experiment we set the rf drive frequency 
at the <5> X stat = nj>n resonant frequency (the pseudoreso- 
nant frequency |5l ll7ll20l |). At this bias flux the ring is at 
its most diamagnetic [with a ring magnetic susceptibility 
X (n&o) = dl s (n§o) /d$ x = — 1] and I s is almost linear 
in $ x . From the data of figure 2 we can make an estimate 
of I c (and hence j3) from the start of the step features in 
the standard SQUID magnetometer compared with the 
(first) step feature in the large scale characteristics poj . 
This yields an I c for the large (3 characteristic = 80uA 
(i.e. a (3 value in the 70 to 80 range). We have shown [24J 
that knowing I c and the Fermi velocity of the material of 
the SQUID ring (here niobium) we can estimate C. Fol- 
ic ~ lOO^A this gives C s ps fewxlO _13 F. In addition, 
from singly connected current-voltage characteristics for 
point contact weak links 19 , a typical value for R s at this 
I c is ps 10^. 

In figure 4 we show the results of solving J5J and 
(J3J to simulate the behaviour of the ring-tank circuit 
system of fig ures 2 and 3, including a 4 Kelvin noise 
source 0, |2(| • Given the very non- linear nature of 
this coupled system, and the number of parameters in- 
volved, an exact match between theory and experiment 
has proved difficult. Nevertheless, we were able to repro- 
duce the important features of the integer bias flux V out 
versus Ii n characteristic of figures 2 and 3 with I c — 80/iA 
and best fit values C 8 = 5 X 10~ 13 F, R s = lOfl and 
f r f — 25MHz. The other circuit parameters are, as for 
figures 2 and 3, Q = 515, K 2 = 0.008, L T = 63nH, 
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FIG. 4: Theoretical (RSJ+C) V out versus 7i n characteristic 
(with a 4 Kelvin tank circuit noise source added) to model the 
data of figures 2 and 3, setting $ xs tat = n& , A s = 6x 10 _10 H, 
it = 6.3 x 10~ 8 H, K 2 = 0.008, Q = 515 and a best fit 
/ r f = 25MHz. Here, the effective triangular ramp amplitude 
modulation frequency is close to 30Hz. 

while the effective triangular ramp frequency (converted 
from an rf current ramp in /^A/sec) is close to 30Hz. As 
can be seen, the computed solution contains both the 
first flat step and the relaxation oscillations on the sec- 
ond (and subsequent) step features of figure 2. Details 
of these (computed) oscillations are shown in the lower 
inset in figure 4. Also shown inset, over a small range 
of Jin in the same figure, is a highly expanded section 
of the first flat step. This too displays (weak) voltage 
oscillations which appear to be either aperiodic or mul- 
tiperiodic. These oscillations were not seen experimen- 
tally, presumably because of limited signal to noise. If 
Ii n is fixed in the regions corresponding to the insets 
in figure 4 - for the experimental data set of figure 2, 
and time series data are collected, these weak oscilla- 
tions can be seen reflected in the (pseudo) phase por- 
trait, as shown in figure 5(a). Here, the pseudo-phase 
portrait is defined as a plot of V ou t (t) against the time 
shifted V ou t (t + t) at all times t, where r is a suitable 
time delay (l0 _4 secin this example). In the main part 
of figure 5(a) we show the pseudo phase portrait, taken 
over a set of complete orbits, for the bias point in the up- 
per inset. In the lower inset we display just one of these 
orbits. It is clear, for the first step, that this pseudo- 
phase portrait simply shows random fluctuations (noise) 
superimposed on a singly periodic orbit. This should be 
contrasted with the pseudo-phase portrait of figure 5(b) 
for Iin biased onto the middle of the relaxation oscilla- 
tions in the second V ou t versus Ii n step feature [shown in 
the upper inset of figure 5(b)]. Here, we have again plot- 
ted a set of orbits, with a single orbit of this set shown in 
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FIG. 5: Experimental pseudo phase portraits (V ou t(t + r), 
taken over sets of complete orbits, versus V (t)) for the two ii n 
bias points of figure 2 as shown arrowed in the insets. Here, 
r = 10 -4 sec and T = 4.2Kelvin. The lower insets in (a) 
and (b) display single orbits of the data sets for the two bias 
points. The data was obtained in a detection bandwidth of 
1MHz with an 18Hz sweep frequency. 



the lower inset. As is apparent both the set and the sin- 
gle orbit display the periodic behaviour associated with 
the relaxation oscillations on the step. 

What the theoretical computations of figure 4, the ex- 
perimental data of figures 2 and 3, and the previously 
published multilevel structures |2£j, demonstrate is that 
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the form of the solutions to the coupled equations of 
motion |J2J and is critically dependent on parame- 
ter values. More specifically, it appears that in the very 
strongly hysteretic regime and above (i.e. j3 > 100) the 
underdamping in the ring-tank circuit system becomes 
sufficiently large to wash out multilevel structures com- 
pletely, or almost completely. Physically, these structures 
arise (with stochastic jumping) from the system moving 
around n hysteresis loops in $ versus "fey (or / s versus 
<& x ) in to rf cycles, n and to integer [22]. When suffi- 
ciently underdamped we assume that the system simply 
jumps between much larger values of $ in <!> versus $ x , 
followed by long recovery period over many rf cycles, i.e. 
a relaxation oscillation process. 

III. CONCLUSIONS 

The experimental data of figures 2 and 3 are typical 
of the V ou t versus Jj„ characteristics found in the very 
strongly hysteretic, and underdamped, regime {f3 > 100). 
These data can be modelled with considerable accuracy 
invoking the SQUID potential and using the quasi- 
classical equations of motion @ and ©, as can be seen 
from the computed V ou t versus Ii n characteristics of fig- 
ure 4. Thus, from experiment and theory we find that 
there exists another class of solution, to the dynamics 
of the driven ring-tank circuit system in the very highly 
hysteretic regime, apart from multilevel structures |20j. 
Since these (latter) structures clearly hold out the real 
possibility of quasi-classical multilevel logic and memory 
if these dynamics are well understood and can be con- 
trolled, there is every reason for opening up the phase 
space of behaviour. Indeed, the potential for single flux 
quantum logic has already attracted a great deal of inter- 
est [3^. It seems reasonable, therefore, that multilevel 
logic based on the manipulation of sets of discrete flux 
states of a SQUID ring in the strongly hysteretic regime 
will also generate interest. From the work reported here, 



it would appear that the change over between multilevel 
and relaxation solutions is quite subtle and sensitive to 
circuit parameters. From experiment and theory we find 
that relaxation oscillation characteristics, very similar to 
those reported here, carry on to very much higher (3 val- 
ues (> 500). This, of course, points to these being a 
rather general class of solutions to J2Jl and J2Jl. Again, 
as with the earlier multilevel structures |20|. these solu- 
tions do not arise in the linearized description of the sys- 
tem. We note, again from experiment and computation, 
that the multi-level structures |20| found in the strongly 
hysteretic regime with (3 ~ 50 gradually fade away as 
(5 increases, to be replaced by the relaxation oscillation 
characteristics reported in this paper. In our opinion, 
even though other factors may be important in govern- 
ing the already underdamped dynamics of this highly 
non-linear system, this points to the damping becoming 
weaker as (3 grows. It also seems clear that a rich new 
field of non-linear dynamics exist which, as yet, has only 
been partially explored. From a historical perspective, 
it is also of interest that the original description of rf- 
biased SQUID magnetometers [2l[ was based on a relax- 
ation oscillation model, albeit operating over very much 
shorter time intervals (a few rf periods), and on the rf 
voltage and current scales of conventional SQUID mag- 
netometers [2]]]. Never the less, relaxation oscillations 
have not been seen in the dynamics of these standard 
magnetometers. To our best knowledge, the relaxation 
oscillations reported in this paper are the first to have 
been observed in the dynamical (V ou t versus 7j„) char- 
acteristics of SQUID ring-resonator systems on any scale 
or frequency. 
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